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Abstract 

' This paper proposes a unified approach to enable the study of diverse distributions 

in the real, complex, quaternion and octonion cases, simultaneously. In particular, 
the central, nonsingular matricvariate and matrix multivariate Pearson type II dis- 
tribution, beta type I distributions and the joint density of the singular values are 
obtained for real normed division algebras. 
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1 Introduction 

In the last twenty years, the concept and the statistical and mathematical techniques known 
as multivariate analysis have changed dramatically. New statistical and mathematical tools 
for the analysis of multivariate data have been developed in diverse areas of knowledge, and 
have promoted new disciplines such as pattern recognition, nonlinear multivariate analysis, 
data mining, manifold learning, generalised multivariate analysis, latent variable analysis 
and shape theory, among others. These and other fields constitute waht is known as modern 
multivariate analysis. 

Renewed interest in multivariate analysis in the complex case has emerged in diverse 
areas, see Metha [26j], Ratnarajah et al. [3(| and Micheas et al. j27|, among many others. 
Similarly, several works involving multivariate analysis have appeared in the context of the 
quaternion case, see Bhavsar [2j, Forrester [15[, Li and Xue [25J, among others. Although 
receiving little attention from a practical standpoint, but equally interesting from the the- 



oretical point of view, some results have appeared in the octonion case, see Forrester 15 1. 
This lack of widespread interest may be, because as asserted by Baez [l[, ...there is still no 
proof that the octonions are useful for understanding the real world. We can only hope that 
eventually this question will be settled one way or another. 
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In terms of concepts, definitions, properties and notation from abstract algebra, it is 
possible to propose a unified approach that enables the simultaneous study of the distribution 
of a random matrix in real, complex, quaternion and octonion cases, which is termed as the 
distribution of a random matrix for real normed division algebras. 

In the real case, the matricvariate Pearson type II distribution appears in the frequen- 
tist approach to normal regression as the distribution of the Studentised error, see Diaz- 
Garci'a and Gutierrez- Jaimez Q and Kotz and Nadarajah 124]. In Bayesian inference, the 
matricvariate Pearson type II distribution is assumed as the sampling distribution; then, 
considering a noninformative prior distribution, the posterior distribution and marginal 
distributions, the posterior mean and generalised maximum likelihood estimators of the pa- 
rameters involved are found, Fang and Li [l3j . A very important question is that of the 
role of the Pearson type II distribution in multivariate analysis, because if the matrix R has 
a matricvariate Pearson type II distribution, then the matrix RR* (or R*R) is distributed 
as beta type I; and the distribution of the latter, in particular, plays a fundamental role in 
the MANOVA model, see Khatri (UHS] and Muirhead [Hj]. 

The present article is organised as follows; a minimal number of concepts and the notation 
of abstract algebra and Jacobians are summarised in Section [5] Section [3] then derives the 
nonsingular central matricvariate Pearson type II and the beta type I distributions and some 
basic properties. Similarly, results are obtained for the matrix multivariate Pearson type II 
and the corresponding beta type I distributions, see Section 2] Finally, the joint densities of 
the singular values are derived in Section [5] We emphasise that all these results are found 
for real normed division algebras. 



2 Preliminary results 

A detailed discussion of real normed division algebras may be found in Baez [l[ and Gross 
and Richards 16!]. For convenience, we shall introduce some notation, although in general 
we adhere to standard notation forms. 

For our purposes, a vector space is always a finite-dimensional module over the field 
of real numbers. An algebra J is a vector space that is equipped with a bilinear map 
m : 5 x 5 5 termed multiplication and a nonzero element 1 £ 5 termed the unit such 
that m(l, a) — m(a, 1) = 1. As usual, we abbreviate m(a, b) = ab as ab. We do not assume 
$ associative. Given an algebra, we freely think of real numbers as elements of this algebra 
via the map lj > uil. 

An algebra J is a division algebra if given a,6 € 5 with ab = 0, then either a = or 
6 = 0. Equivalently, 5 is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra J that is 
also a normed vector space with ||ab|| = ||a||||6||. This implies that J is a division algebra 
and that ||1|| = 1. 

There are exactly four normed division algebras: real numbers (5ft), complex numbers 
(£), quaternions (Sj) and octonions (£)), see Baez We take into account that 5ft, £, f) and 
D are the only normed division algebras; moreover, they are the only alternative division 
algebras, and all division algebras have a real dimension of 1, 2, 4 or 8, which is denoted by 
/?, see Baez [l|, Theorems 1, 2 and 3]. In other branches of mathematics, the parameters 



a = 2/ j3 and t — /3/4 are used, see Edelman and Rao [12| and Kabe [20(, respectively. 

Let £^ „ be the linear space of all m x n matrices of rank m < n over 3" with m distinct 

positive singular values, where $ denotes a real finite- dimensional normed division algebra. 

Let g rmx " be the set of all m x n matrices over 5- The dimension of fi mxn over 5ft is ftmn. 

t 

Let A £ ffnxn^ thgjj a* = A denotes the usual conjugate transpose. 

Table [T] sets out the equivalence between the same concepts in the four normed division 
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algebras. 



Tabic 1: Notation 



Real 



Complex 



Quaternion 



Octonion 



Semi-orthogonal Semi-unitary Semi-symplectic 



Orthogonal 
Symmetric 



Unitary 
Hermitian 



Symplcctic 

Quaternion 
hermitian 



Semi-exceptional 
type 
Exceptional 
type 
Octonion 
hermitian 



Generic 
notation 
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In addition, let be the cone of positive definite matrices S € y™ xm ; then is an 
open subset of 6^. 

Let 2)^ be the diagonal subgroup of Cf^ m consisting of all D 6 ^ mxm ; D = diag(di, . . . , d m ). 

For any matrix X £ ^nxm^ ci e xiotes the matrix of differentials (dxij). Finally, we 
define the measure or volume element (dX) when X £ S mXTl , 6^, or n , sec Dimitriu 
Hi- 

If X € ynxn ^ nen ((jX) (the Lebesgue measure in g ,mxn ) denotes the exterior product 
of the fjmn functionally independent variables 



(dX) = A A dxy where dajy = A dx 
i=i i=i 



(fe) 



fc=i 



If S e 6^ (or S e %{{m) is a lower triangular matrix) then (dS) (the Lebesgue mea- 
sure in or in T^(m)) denotes the exterior product of the m{m + l)/3/2 functionally 
independent variables (or denotes the exterior product of the m(m — l)/3/2 + n functionally 
independent variables, if su € 5ft for ell i = 1, . . . , m) 



,(*) 



A A *l 

(dS) = <( ''-^ fc = 1 o 

/\ds u /\ /\*»«\ ifs «e5R. 

i=l i<J fe=l 

The context generally establishes the conditions on the elements of S, that is, if £ 5ft, 
£ £, £ ij or € D. It is considered that 

ra {3 m m 

( dS ) = A A*# )s A*«A A**- 

i<j fc=l i=l i<j fe=l 

Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S € ^3^, 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive matrix). 

If A £ 2)^ then (dA) (the Legesgue measure in 23^) denotes the exterior product of the 
f3m functionally independent variables 



n 



(dA) = A A dA 



i=i k=i 



If Hi e _ then 



(H*dHi) 



A A h ;^- 
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where H = (H*|H?,)* = (hi, . . . , h m |h„ l+ i, . . . , h„)* G il^(n). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When n = 1; j defines 
the unit sphere in $ m . This is, of course, an (m — l)/3- dimensional surface in $ m . When 
n = m and denoting Hi by H, (HdH*) is termed the Haar measure on iiP(rn). 
The surface area or volume of the Stiefel manifold n is 

Vol«J = / (HidHt) = (1) 



where r^Ja] denotes the multivariate Gamma function for the space 6^, and is defined by 



/ etrl-AHAI -^^ 1 ^/ 2 ^ 1 ^) 

m 



i=l 



where etr(-) = exp(tr(-)), | • | denotes the determinant and Re(a) > (m — l)/3/2, see Gross 
and Richards flq . Similarly, from Herz [l8j the multivariate beta function for the space 
6^, can be defined as 

B&M = / | B | a -c- m -i)^/2-i| Im _ B | b -( m+1 ^/ 2 - 1 (dB) 

Jo<B<I m 

r&[«]r&[6] 



r&[a + 6] 



(2) 



where A = (I - B)- 1 - I, Re(a) > (m - 1)0/2 and Re(6) > (m - l)/3/2. 

Now, we show three Jacobians in terms of the /3 parameter, which are based on the work 
of Kabe and Dimitriu These results are proposed as extensions of real, complex or 
quaternion cases, see James |l9j|. Khatri (22[, Metha {2(1, Ratnarajah et al. [3(j and Li and 
Xue 25j, also see Diaz-Garcia and Gutierrez- Jaimez [6|. 

Lemma 2.1. Let X and Y g C^ n , and let Y = AXB + C, where A e £^ m , B G £^.„ 
and C € £^ n are constant matrices. Then 

(aY) = |A*A|' 3n / 2 |B*B| /3m / 2 (dX). (3) 

Lemma 2.2 (Singular value decomposition, SVD). Let X £ Cf^ n , such that X = V'DWi 
with V G U^m), Wi G „ and D = diag(di, • • • , d m ) € S^, di > • • • > d m > 0. TTiera 

m m 

(dX) = 2" m ^ T ]^[df ( ™" m+1K1 [|(d 2 -d 2 ) /3 (dD)(VdV*)(WidWJ), (4) 

<=i 

w/iere 







< o, 


= 1; 


J "TO, 


= 2; 


] ~2m, 


= 4; 


[ -4m, 


= 8. 


XX* G 


g,. TTien 



Lemma 2.3. Let X G £^ n , and ! 

(dX) = 2- m |S|^ ,l - m+1 )/ 2 ~ 1 (dS)(V 1 dVJ). (5) 
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3 Matricvariate Pearson type II distribution 

In the real case, the matricvariate Pearson type II distribution (also known in the literature 
as matricvariate inverted T distribution) was studied in detail by Dickey II Oil and 
Cadet [3|, see also Press [29(. This distribution was previously studied by Khatri [21[, also 
in the real case. 

Theorem 3.1. Let R e n defined as 

R= L _1 X 

where L is any square root of U = (Ui + XX*) such that LL* = U, Ui ~ W^(f, H), 
v > /3(m-l), independent o/X - A/^ xn (0, l m ®B), and H € T/ien U - W&(i/+ra, If) 
independently of R. Furthermore, the density of R is 

rgJ fi n tr^^l I 1 ™ ~ HRI*"-"^- 1 , I m -RR*e^ (6) 

which is termed the matricvariate Pearson type II distribution^. 

Proof. From Kabe [2(| and Diaz-Garcia and Gutierrez- Jaimez 6, 9] the joint density of Ui 
and X is 

oc lUil^-^^/^et^-jSH-^i+XX*)^}, 
where the constant of proportionality is 

1 1 



(2/3-i)^/2r^[/3i//2]|H|^/ 2 (2 7 r/3- 1 )' 3mn /2|H|' 3 «/ 2 ' 

Making the change of variable Ui = (U — XX*) and X = LR, where U = LL*. then by {3} 

(<*Ui)(dX) = | LL* | ^"Z 2 (rfU) (dR) = |U|^/ 2 (dU)(dR), 

and observing that |Ui| = |U-XX*| = |U-LRR*L*| = |U||I m -RR*|, the joint density 
of U and R is 

oc |I m - RIi *^- m +i)/2-^ v ^+n- m +i)/2-i etr{-^H- 1 U/2}. □ 

□ 

Similarly to Dickey 10], ([6]) is alternatively given by 
Ti\p{n + v)/2] 



n mn P/ 2 rP[/3(n + v-m)/2] 



|I„ - R*R|fl("- TO + 1 )/ a - 1 ) I„ - R*R e (7) 



Corollary 3.1. Let Q = (M*) _1 RN _1 + /x, R as in Theorem \3.1l M and N are any 

square root of the constant matrices B— MM* e and JJ= NN* € *p^, respectively, and 
H G £^ n is constant. Then, from {6]) and f?p £/ie density of Q is given by 



1 In the literature it is customary to use the real matricvariate Pearson type II distribution, complex 
matricvariate Pearson type II distribution, quaternion matricvariate Pearson type II distribution and octo- 
nion matricvariate Pearson type II distribution,; here, however, we use simply matricvariate Pearson type 
II distribution as the generic term. 
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fl- i -(Q- / i) 1 zr 1 (Q-Mre?p& > 

and as 

^[P{n + v)/2] \B\M , (v _ m+1)/2 _ 1 

^^(,| n - m )/2]W^ 1)/2 - lW ^ M) ^ M)l 

^-(q-m)*b(Q-m) e?pg. 

T7us /aci shall denoted as Q ~ VIIj nxn {v, /x, 5, ,ZZj fane! 0/ course R ~ VII^ lxn (i / , 0, I m , 
Proof. The proof follows observing that, by (|3|) 

(dR) = \MM*f n/2 \NN*f m/2 (dQ) = |B|^"/ 2 |/i;| /3m/2 (dQ) 

and that 

IB- 1 - (Q - m)A _1 (Q - = IBI^IAI^IA - (Q - M)*B(Q - fi)\. □ 

□ 

Corollary 3.2. Assume that Q ~ "P//f lx „(V, /x, B,/f), then 

Q* ~ VII% xm (v, fi* , B' 1 )- 

Proof. The proof follows immediately from the two expressions for the density of Q in 
Corollary O □ □ 



Corollary 3.3. Let R € C^ n n defined as 

R = YL^ 1 

where Li is any square root of V = (Vi + Y*Y) smc/i i/iai LiLJ = V. Vi ~ (y + n — 
m, .77), independent of Y ~ A/^ x „(0, J® I„), and J6 ^. Traen V - W£0 + ra, .ZZ) 
independently o/R. Furthermore, R ~ 'P/7^ xn (i/, 0, J m , J n ). 

Proof. The proof is a verbatim copy of the proof of Theorem 13.11 □ □ 

Now, assume that R ~ 7 , //^ xn (f, 0, I m , I„) with n > m and let B g defined as 
B = RR* then, under the conditions of Theorem 13.11 and Corollary 13. 31 we have 

B = L _1 XX*(L -1 )* = L- 1 W(L" 1 )* 
= Y(Vi + Y*Y)~ 1 Y*, 

where W = XX* - W&(n, H), n > /3(m - 1). Thus: 
Theorem 3.2. The density o/B is 

1 | B |^(n-m+l)/2-l| I _ B |/J( v _m+l)/2-l /g) 

^[^/2,/3n/2]' 

where B^ n [-, ■] is given by 0) and B is said to have a matricvariate beta type I distribution. 
Proof. The proof follows from ([6]) by applying (JlJ and ([5]). □ □ 
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In addition, assume that n < m and let B £ denned as B = R*R then, under the 
conditions of Theorem 13.11 and Corollary 13. 31 we have 



B = X*(Ui +X*X)~ 1 X 

where W, = Y*Y - W^(m,H), m> fi(n- 1), Thus: 
Theorem 3.3. B has the density 

1 



BZ\J3(v + n-m)/2,Pm/2] 
Also, we say that B has a matricvariate beta type I distribution. 



B |/3(m-™+l)/2-l|j n _ B |/3( I /-m+l)/2-l_ ( g ) 



Proof. The proof is the same as that given in Theorem l3.2l Alternatively, observe that den- 
sity © can be obtained from density ([8]) making the following substitutions, see Muirhead 



28|, Eq. (7), p. 455] and Srivastava & Khatri [3lJ, p. 96], 

m — > n, n — > m, v v + n — m. □ (10) 

□ 

Densities (jHJ and (j9]) have been studied by several authors in the real case, see Khatri 
23] and Srivastava & Khatri (3l|, Cadet Q, among many others; and by James [3], Dfaz- 
Garci'a and Gutierrez-Jaimez [7j and Diaz-Garcia and Gutierrez-Jaimez [5j and Di'az-Garci'a 
and Gutierrez-Jaimez Q, in noncentral, doubly noncentral, singular and nonsingular and 
complex cases, among many other authors. By the analogous construction of B = RR* (or 
B = R*R) in terms of R compared with the construction of the Wishart matrix in terms 
of matrix multivariate normal distribution, the distribution of B = RR* (or B = R*R) is 
sometimes termed the studentised Wishart distribution. Kabc 20] studied densities ([8]) and 
© for the hypercomplex case. 

4 Matrix multivariate Pearson type II distribution 

Theorem 4.1. Let (5 1/2 ) 2 ~ X 2 ^(v) independent of Y ~ J\ff nxn (0,I m (g> I n ), and de- 
fine Si = S + tr YY* and Ri = S t 1//2 Y. Then Si <~ x 2 (^ + ran) independently of Ri. 
Furthermore, the density of Ri is 

r^(, + mn)/2] M2 _ 1-trlWUX), (11) 

which is termed the matrix multivariate Pearson type II distribution. 
Proof. The joint density of S and Y is 



CX S' 



\uj1-X 



etr{-/3(s + trYY*)/2}, 



and the desired results are obtained analogously to the proof of Theorem 13.11 □ □ 

Corollary 4.1. Let Qi = (M*) _1 RiN _1 + fi, Ri as in Theorem HOI M and N are any 

square root of the constant matrices B— MM* £ and NN* £ respectively, and 
fi £ n is constant. Then, 

lfX + r ^!m W»IM'"'» [1 - tr E(Q, - „yBQ, - rtf" . 
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1 - tr E(Q 1 - n)MQi - M)* >0- 

Hence, we write 

Qx~MPII^Jy,n,B,A), 
and, in particular R4 ~ M.VH^ nxn (v, 0, I m , I„). 
Proof. The proof follows observing that, by ((3|) 

(dRi) = |MM*|^" /2 |NN*| /3m/2 (dQ 1 ) = |B|^ n/2 |A|^ m/2 ((iQi)- □ 

□ 

Now, assuming that Ri ~ VII^ lxn (y, 0, I m , I n ), with n> m and defining Bi = RiR* € 
^P^j, then, under the conditions of Theorem 14. II we have that 

Bi = S^YY* 

where W = YY* ~ W^(n, I m ), n > f3(m — 1). Furthermore: 
Theorem 4.2. The density of Hi is 



rjg^+mn)/2] ^ (n _ m+1)/2 _ 1 _ )M 2-i B j 
rf[#//2]lt[/3n/2] 



(12) 



Bi is said to /iave a matrix multivariate beta type I distribution. 

Proof. The proof follows from (fTTj) by applying (TTJ) and (|5) . □ □ 

Similarly, if n < m and Bi = R|Ri 6 *p£, thus: 
Theorem 4.3. Bi ftas i/ie density 



Y{\fi(y + mn)/2\ ~ |/3(ro _ n+1)/2 _ 1(1 _ | ^/j-i „ g j 
lf[^i//2]r^m/2] 



(13) 



Thus, Bi is sairf to Ziaue a matrix multivariate distribution type I distribution. 
Proof. Density (fT3")) can be obtained from density ([12")) by making the following substitutions, 

m — > n, n — > m. □ (14) 

□ 

5 Singular value densities 

In this section, we derive the joint density of the singular values of matrices R, R, Ri and 
Ri. In addition, the joint densities of the eigenvalues of B, B, Bi and Bi are obtained. 

Theorem 5.1. Let Si, . . . , S m be the singular values of R ~ T > II^ nxn {i', 0, I m , I„), 1 > 5\ > 
■ ■ ■ > 5 m > 0. Then its joint density is 



T p m [f3m/2}B^/2,f3n/2}l\ , . 



—————— n s^- m - d»)«-™+D/»-i - (15) 



where r is defined in Lemma 



S 



Proof. The proof follows immediately from © , first using Q and then applying (TI| . □ □ 

The joint density of the singular values of R is obtained from (jT5j) after making the 
substitutions (flU)) . 

Theorem 5.2. Suppose that Ri ~ AdVIIm Xn (i>, 0, I m , I„) and Zei ai,...,a m , 1 > a\ > 
■ ■ ■ >,a m > 0, < J2iLi a i < 1; Ms singular values. Then its joint density is 



Tf^/2]T^m/2]T^\pn/2] 



r^r^./o^r^/o^r^/oi 1 ~ H a * lK«i " ( 16 ) 



i— 1 / z— 1 i<j 



Proof. The proof is analogous to that given for Theorem 15.11 □ □ 

Similarly, the joint density of the singular values of Ri is obtained from (|16p and making 
the substitutions ([H]). 

Finally, observe that Si — ^/eig^RR*) and en — ^/eig^RiRJ), where eig,(A), i = 
1, . . . , m, denotes the i-th eigenvalue of A. Let Ai = eigj(RR*) and % — eig,(RiR|), hence, 
observing that, for example, Si — y/Xl and then 

m mm 

/\dSi = /\2- m Y[x; 1/2 d\ i , 

i—l i—1 i—l 

the corresponding joint density of Ai, . . . , A TO , 1 > Ai > • • • > A m > is obtained from ([T5]l 
as 



T m [pm,/2}Bm{f3v/2,f3n/2} 



J] Af ("- m + 1 )/ 2 - 1 (l - A .)/^-m+l)/2-l JJ (A . _ Aj .y 



Tfl™ /olH/3 

■ II// ' I i I '■' . " I 



Similarly, the joint density of 71, . . . ,7 m , 1 > 71 > • • • > J m > 0, < YnLi a i < 1> i s 
obtained from (fT6|) as 

Conclusions 

• Beyond a doubt, in any generalisation of results there is a price to be paid, and in this 
case the price is that of acquiring a basic understanding of some concepts of abstract 
algebra, which can initially be summarised as the use of notation and a basic minimum 
set of definitions. However, we believe that a detailed study of mathematical properties 
from a statistical standpoint can have a potential impact on statistical theory. 

• Furthermore, note that X £ £^ n has a matrix multivariate elliptically contoured 
distribution for real normed division algebras if its density, with respect to the Lebesgue 
measure, is given by (see Di'az-Garci'a and Gutierrez- Jaimez jy]): 

WM^l/. h ^ [^( X - ^ X ( X - /*)*] } ■ ( 17 ) 

where n £ £f n „, S £ <£f n , © £ <£f n . The function h : $ -> [0, 00) is termed the 
generator function, and it is such that f^p u l3nm ~ 1 h(u 2 )du < 00 and 

cv,»)=^M{/ ¥ y""-'% 2 >4 
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Such a distribution is denoted by X ~ £„. rm (u, S, 0, h), for the real case, see Fang 



and Zhang [1J] and Gupta, and Varga [ljj; and Micheas et al. |27[ for the complex 



case. Observe that this class of matrix multivariate distributions includes normal, con- 
taminated normal, Pearson type II and VII, Kotz, Jensen-Logistic, power exponential 
and Bessel distributions, among others; these distributions have tails that are more or 
less weighted, and/or present a greater or smaller degree of kurtosis than the normal 
distribution. 

Assume that X = ( Xi : X2 ) ~ £^ ixn+v (0,l m ,l n + 1/ ,h), n,v > m; and define, 

rn X n rnxu 

R = L _1 Xi, where L is any square root of V = (X 2 X 2 + XiX*) such that LL* = V. 
Then R - VII^ xn (v,0,I m ,ln) independently of V - QW^n + u,I m ,h), n + v > 
/3(m — 1), where SW^(-) denotes the generalised Wishart distribution based on an 
elliptical distribution, see Dfaz-Garcfa and Gutierrez- Jaimez Q and Diaz-Garcia and 
Gutierrez- Jaimez Q . From (jTTJ) the density of X is 

CP(m, n + v)h {tr (XiXJ + X 2 X*)} . 

Let V = X 2 X*, then by ([23]) , (dX 2 ) = 2- m |V | /3(l/ - m+1)/2 - 1 ( dV o)(HidH^. Thus, 
the marginal density of Xi and Vo is obtained by integrating over Hi £ „, and 
then using ([1]), to obtain 

C ' K :r!''"" |v.l"-""»>l«. (x,x; + v„» . 

Now, let V = (V + XiXJ) and R = L- : Xi, where LL* = V, then by (|2T|) 

(dX x )(dVo) = \Vf n / 2 (dR)(dV). 
Observing that | Vol = |V||I m — RR*|, and so the joint density of R and V is 

- RRf M+')/M|yf (^-n+WM,,{ tl V}, 

from where the desired result follows. □ 

Observe that in this case, Xi and X 2 (or Vo = X 2 X 2 ) are stochastically dependent. 
Furthermore, note that only when the particular matrix multivariate elliptical distri- 
bution is the matrix multivariate normal distribution, are Xi and X 2 (or Vo = X 2 X 2 ) 
independent. Then, we can say that the matricvariate Pearson type II distribution is 
invariant under the family of matrix multivariate elliptical distributions, and further- 
more, its density is the same as when normality is assumed. In the same way, it can be 
proved that the matrix multivariate Pearson type II, matricvariate and matrix multi- 
variate beta type I distributions are invariant under the family of matrix multivariate 
elliptical distributions. 

• Finally, following Kabe [2(| , the distributions studied in this paper are easily extended 
to the hypercomplex case (biquaternion and bioctonion cases, which are a Jordan 
algebras), simply replacing j5 by 2/3 in the obtained results. 
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